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Hadronic decays from the lattice 

UKQCD Collaboration: C. Michael 8 and C. McNeile a 

a Theoretical Physics Division, Dept. of Mathematical Sciences, 
University of Liverpool, Liverpool L69 3BX, UK 

We discuss strategies to determine hadronic decay couplings from lattice studies. As an application, we explore 
the decay of a vector meson to two pseudoscalar mesons with Nf = 2 flavours of sea quark. Although we are 
working with quark masses that do not allow a physical decay, we show how the transition rate can be evaluated 
from the amplitude for p — > nn and from the annihilation component of -kit — > nn. We explore the decay amplitude 
for two different pion momenta and find consistent results. The coupling strength we find is in agreement with 
experiment. We also find evidence for a shift in the p mass caused by mixing with two pion states. 



1. INTRODUCTION 

The recent renaissance of experiments to study 
the spectroscopy of light hadrons is partly driven 
by the goal to understand confinement. Any anal- 
ysis of experimental data requires both a calcu- 
lation of decay widths as well as the masses pp. 
Although computing the masses of resonances 
is (currently) not part of the program to com- 
pute "gold plated" observables with high preci- 
sion from lattice QCD, dealing with decay widths 
is an inherent part of studying hadrons which are 
under current experimental scrutiny, such as "ex- 
otics" , and gluonic components of scalar mesons. 
Here we compute the coupling of the p meson to 
two pions to validate our procedure. 

One limitation of the lattice approach to QCD 
is in exploring hadronic decays because the lat- 
tice, using Euclidean time, has important con- 
tributions from low lying thresholds which 
can obstruct the study of decay widths. The fi- 
nite spatial size of the lattice implies that two- 
body states are actually discrete. By measuring 
their energy very precisely as the spatial volume 
is varied, it is possible to extract the scat- 
tering phase shifts and hence decay properties. 
For on-shell transitions, it is possible to estimate 
hadronic transition strengths more directly and 
this approach has been used to explore 0] hybrid 
meson decay rates. Here we explore this approach 
further for the case of p meson decay to 7T7t, fol- 
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2. TRANSITIONS 

The situation we shall analyse is represented by 
the energy spectrum shown in fig. Q here neglect- 
ing any interactions among the states. We eval- 
uate correlations between lattice operators cre- 
ating both a p meson and a 7T7t state, using a 
stochastic method to evaluate the quark diagrams 
shown in fig. [5] from 20 gauge configurations 6 
(double the number used in our initial calcula- 
tion) with Nf = 2 sea quarks of mass corre- 
sponding to about 2/3 of the strange-quark mass. 
These correlations, normalised by the two point 
functions as appropriate, are illustrated in fig. |3| 
Here the off-diagonal case (labelled p\ — > ttittq) 
shows the important feature that it grows approx- 
imately linearly with increasing t. This linear 
growth will only occur for on-shell transitions 4 
and this is essentially the case here. 

To extract an estimate of the transition am- 
plitude, we introduce a parameter x = {p\tttt) 
where these states are normalised on the lattice 
(to unity) . If higher excited states are neglected, 
one can make a two state model (with basis states 
p and the lightest 7T7t state) with this transition 
amplitude and evaluate the contributions shown 
in fig. Unlike weak decays, there is no specific 
operator that causes the transition in strong de- 
cays, hence the specific time (between and t) of 
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Figure 2. 



any transition, such as p — * 7T7t, is not known. In 
general this makes it hard to control the contri- 
butions from excited states. As emphasized pre- 
viously these excited state contributions can 
be avoided if the transition is approximately on- 
shell, when the linear dependence on t of the off- 
diagonal transition amplitude (with the slope x) 
is a unique signature of this on-shell transition. 

Note that for the case of po — > ttitti, the on- 
shell condition is much less well satisfied but the 
relative momentum of the pions in the centre of 
mass is twice as large as for pi — > n ni and hence 
x should be approximately twice as large, since 
for a P-wave decay there will be a momentum 
factor in the transition amplitude. 

A further check of the extraction of x comes 
from the box diagram, fig. which will have a 
contribution behaving as x 2 t 2 /2 arising from a p 



intermediate state - see fig. |3] 
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Figure 3. Ratios of 3 and 4-point correlators. 



For a more quantitative estimate of x, one can 
suppress excited state contributions in analysing 
the data of fig. El as discussed in ref. p] . 

Note that the above methods could be used, in 
principle, even in the quenched approximation. 
They depend on assuming that xa is fairly weak, 
as we indeed find. 

3. ENERGY SHIFTS 

A more rigorous approach is to focus on energy 
values. When two levels are close (our on-shell 
condition), then they will mix and the resultant 
energy shifts give relevant information [3j . More- 
over we can estimate these energy shifts from 
our x- value which provides more cross checks. 
These shifts can only be studied using dynami- 
cal fermions. 

From a full variational analysis we obtain the 
energy shift of the ttittq state (i.e. the un-binding 
energy), as needed in Luscher's approach, as 
0.02(2) upward which is consistent but not suffi- 
ciently accurate to use. The energy shift of the pi 
state can however be determined because of a lat- 
tice artifact. The p with momentum 1 (in lattice 
units of 2n / L) can have its spin aligned parallel 
to the momentum axis (P) or perpendicular to it 
(A). Because the ttqtti state has relative momen- 
tum along a lattice axis and the transition from p 
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to 7T7T has orbital angular momentum L=l (so a 
distribution like cos 8), only the parallel state (P) 
can mix with this two pion state. This mixing will 
not be present in the quenched approximation, so 
this provides a direct opportunity to see the effect 
of the two pion channel on the p in unquenched 
studies. We do indeed find such a mass splitting 
between the P and A orientations of the p due to 
mixing for Nf = 2 but not for Nf — as shown 
in fig. 0] Moreover the magnitude of this energy 
shift (0.026 ±7 in lattice units) is consistent with 
other determinations of the transition strength. 
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Figure 4. The ratio of Parallel to Antiparallel p 
correlators. 

4. PHENOMENOLOGY 

The basic assumption is that the transition 
from p to 7T7T is given by an effective interaction 
with a finite spatial extent, this is usually sum- 
marised by an effective lagrangian where we nor- 
malise the coupling as g 2 — TME/k 3 in terms of 
the decay width. Then, provided the lattice spa- 
tial size is big enough that the hadrons are not 
distorted, our lattice situation (where no decay 
occurs) can be used to determine g and this can 
then be used to predict decay widths when quark 
masses are varied, assuming that the coupling is 
largely independent of the quark masses. 

From our lattice studies, we deduce that ax — 
0.06^ for pi —* ttittq. Translating |S] this lattice 



transition amplitude to the continuum normali- 
sation, gives g = 1.40^23 • Using the observed 
pi energy shift gives another estimate, namely 
g = 1.56^3- Note that our lattice values would 
need to be extrapolated to light sea-quarks and to 
the continuum limit to allow all sources of system- 
atic error to be explored. Nevertheless, these two 
values agree well with the values extracted from 
decays of p, K* and 4> mesons, namely g « 1.5. 

Encouraged by the success of the calculation 
of the decay width of the p meson, we are start- 
ing to look at the decays of scalar mesons on the 
lattice. This is important for studying the glue 
components of singlet scalar mesons and to study 
ao (980) and /q (980) mesons which are believed to 
have important 4 quark contributions. There has 
been one earlier calculation of the decay width of 
the ++ glueball from lattice QCD, but this has 
not been exploited by phenomenology yQ. 

For decays of singlet scalar mesons there are ad- 
ditional diagrams to those in figure[5]from discon- 
nected graphs. The "noise" from the additional 
diagrams makes the extraction of a signal hard. 
In the real world the non-singlet ao(980) meson 
has the strong decays ao — > ttt] and an — > KK. A 
direct comparison with experiment will be non- 
trivial because the nrj decay involves a noisy dis- 
connected loop and the ao — > KK decay requires 
Nf = 2 + 1 sea quarks for a theoretically clean 
calculation in this channel. 
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